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Abstract 

Results of FLMTO-GGA (full-potential linear muffin-tin orbital - generalized gra- 
dient approximation) calculations of the band structure and boron electric field gra- 
dients (EFG) for the new medium-T c superconductor (MTSC), MgB2, and related 
diborides MB2, M=Be, Al, Sc, Ti, V, Cr, Mo and Ta are reported. The boron EFG 
variations are found to be related to specific features of their band structure and 
particularly to the M-B hybridization. The strong charge anisotropy at the B site in 
MgB2 is completely defined by the valence electrons - a property which sets MgB2 
apart from other diborides. The boron EFG in MgB2 is weakly dependent of ap- 
plied pressure: the B p electron anisotropy increases with pressure, but it is partly 
compensated by the increase of core charge assymetry. The concentration of holes 
in bonding a bands is found to decrease slightly from 0.067 to 0.062 holes/B under 
a pressure of 10 GPa. Despite a small decrease of N(Ejd, the Hopfield parameter 
increases with pressure and we believe that the main reason for the reduction under 
pressure of the superconducting transition temperature, T c , is the strong pressure 
dependence of phonon frequencies, which is sufficient to compensate the electronic 
effects. 
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Recently, Akimitsu et all reported the discovery of medium-T c superconductivity 
(MTSC) with T c of about 39 K in magnesium diboride (MgB 2 ) with a simple composi- 
tion and crystal structure (AlB2-type, space group P6/mmm, Z=l). The band structure 
calculations showedEHl that the MTSC in MgB2 can be attributed to a strong electron- 
phonon coupling, a rather high density of states from 2D (in-plane) metallic boron cr(p Xty ) 
bands at and the existence of p^-band holes. By now, a number of studies have been 
performed with NMrJHS which is a very powerful technique for investigating the properties 
of MgB 2 . The measured quadrupole interaction is determined by the value of the electric 
field gradient (EFG), which is directly related to the charge distribution around the nucleus. 
Thus, theoretical EFG studies are important in order to give a reliable interpretation of the 
experimental data based on the electronic structure. 

The EFG at the boron site in MgB 2 was found experimentally to be much larger than 
those for d- diborides (Table I). As seen from Table I, the variation in EFG for diborides 
covers two orders of magnitude and shows trends which cannot be explained by the crystal 
structure changes. For example, these EFG's in A1B 2 and TiB 2 differ by almost three times, 
but their lattice parameters a and c are approximately the same; and, vice versa, equal 
EFG's were obtained in BeB 2 and MgB 2 , for which the lattice parameters have the largest 
differences among diborides under consideration. No attempts was made previosuly to relate 
the EFG changes to the peculiarities of electronic structure. 

In this paper, we present results of first-principles full-potential LMTO- GGA0 calcula- 
tions (within the generalized gradient approximation (GGA) for the exchange correlation 
potential) of the electronic structure and EFG's at the boron site for MgB 2 and other s-,p- 
,<i-diborides. We compare the calculated EFG's with the experimental and other theoretical 
data and explain the EFG's variation based on the anisotropy of boron 2p occupancies. 
The pressure dependence of the EFG, which is very sensitive to the charge distribution, 
represents a good test for the anisotropy study. The experimental and theoretical data on 
the electronic and elastic behavior of MgB 2 under pressure are contradictory: it was found 
to be nearly isotropicSE!, anisotropic^ or strong anisotropic^^. We simulated the pres- 



sure effect on the band structure of MgB2 and estimated the changes in the EFG, boron 
p-occupancies, hole concentration and Hopfield parameter under pressure. These investiga- 
tions are of interest, since the pressure dependence of the T c is the key difference between 
the conventional BCsi and hole0 superconductivity mechanisms. 

The band structures of MgB 2 and some other diborides are shown in Fig.l. There 
are two distinct sets of B 2p-bands: a (2p XjV ) and it (j» 2 )-types with considerably different 
dispersions. The B 2p XjV bands are quasi-two dimensional (2D) along the T-A line in the 
Brillouin zone (BZ) and make a considerable contribution to the density of states at the 
Fermi level, N(E,p), for MgB2. Now was shownUl, that the existence of degenerate p x ,y- 
states above Ep at the T point in the BZ is crucial for the MTSC in diborides. The high Tc 
is explained by the strong coupling of these holes to the in-plane E 2ff phonon mo deffl. The 
B 2p z -bands are responsible for the weaker pp^-interactions and these 3D-like bands have 
maximum dispersion along T — A. The bonding and antibonding B p z bands cross E^ at 
the K point and their location and dispersion depend on the M-B hybridization. For BeB 2 
(Fig. lb) and A1B 2 (Fig.lc), the p XtV bands are, respectively, partly and completely filled, 
the Fermi surface topology changes and medium -T c superconductivity is absentffl. The 
p z -bands progressively move down in going from BeB 2 to MgB 2 and A1B 2 , demonstrating 
the strengthening of M-B bonding. 

The band structure and chemical bonding of all 3d, Ad and 5d-metal diborides were 
previously investigated in detailf^HEI. These studies showed that the cohesive properties of 
AlB 2 -type diborides are explained in terms of the band filling. The Fermi level for TiB 2 (ZrB 2 
and HfB 2 ) falls in the pseudogap where bonding states are occupied and antibonding states 
are empty (Fig. Id). The shift of E^ results in the partial emptying of the bonding states 
(ScB 2 , Fig.le) or the occupation of antibonding states (VB 2 , CrB 2 , MoB 2 , TaB 2 ). Both 
cases correspond to the lowering of the cohesive properties (melting temperature, enthalpies 
of formation, etc.). From first-principles estimates of M-M, M-B and B-B bonding strengths, 
we found that the cohesive energy of <i-diborides with filled bonding states decreases when 
the atomic number increases in the row due to the weakening of M-B hybridization and 



that the Allotype diborides of the group VII and VIII elements are unstableEHi. Thus, 
the main feature of the band structure of TM diborides is the progressive filling of B p x ,yi 
p z -bands and the dominant role of d-states near E^. 

A systematic search for superconductivity in the <i-diborides (with M = Ti, Zr, Hf, V, 
Ta, Cr, Mo) showed that T c is below ~ 0.4 K0 Only NbB 2 and ScB 2 were found to be 
superconductors with a T c of about 0.6 K0 and 1.5 K@, respectively. Recently, a relatively 
high critical temperature, T c ~ 9 K, was found in TaB 2 0. As seen in Fig. If, the bonding 
states in TaB 2 are fully occupied, the Fermi level is shifted away from the pseudogap to 
the region of antibonding states and Ta 5d states define N(E^) (0.9 states/eV). Among the 
3<i-diborides, only in ScB 2 are the 2p x ^ y bands not completely filled (there is a small hole 
concentration of these states at A), but they lie below E^ at T and the largest contribution to 
N(E^) arises from Sc 3d states. Based on the band structure results and calculations of the 
electron-phonon interactional, one may conclude that superconductivity with medium T c is 
unlikely in undoped diborides, except for MgB 2 ; the absence of holes in the two-dimensional 
cx-bands at T results in hardening of zone-phonon modes and weak electron-phonon coupling. 

The electric field gradient tensor, defined as the second derivative of the electrostatic 
potential at the nucleus, was calculated directly from the FLMTO charge density. The 
calculated principal components of the boron EFG tensor, V£, are shown together with 
other theoretical and experimental data in Table 1. Note that the asymmetry parameter, 
(I Vxx I — l^j/j/ D/I^zzl; is equal to for the AlB 2 -type structure. 

The largest boron EFG's in MgB 2 and BeB 2 demonstrate the strongest assymetry of the 
charge distribution as compared with other diborides (note, that here and below we consider 
the absolute value of EFG's). For the 3<i-diborides, the EFG decreases from ScB 2 to TiB 2 
and increases when going from TiB 2 to VB 2 and CrB 2 . The boron EFG's for Ad (MoB 2 ) and 
5d (TaB 2 ) diborides are much smaller than the EFG's for isoelectronic 3d diborides. Note 
that all calculated EFG's are in very good agreement with available experimental data and 
with FLAPwil theoretical results!! (Table 1). 

To analyze the variation of the boron EFG's in diborides, we consider it as a sum of 
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electron (V zz ) and lattice (ion) (V z f) contributions (Table 1). For the s,p diborides, the 
ion contribution is small, and the boron EFG is mainly determined by the anisotropy of the 
valence electrons. For the d— diborides, V\f is larger and depends somewhat on the metal, 
(except ScB 2 where the lattice parameters are the largest). For TaB 2 , the electron and ion 
contributions are almost equal, the EFG is positive and smallest among all diborides. The 
electronic contributions explain the boron EFG variation, although they overestimate the 
calculated and observed EFG's; clearly, the lattice contributions must be taken into account 
in order to obtain good agreement with experiment for d diborides. Among the diborides, 
the lattice contribution is the smallest for MgB 2 and the strong charge anisotropy in MgB 2 
is completely defined by the valence B p electrons - a property which sets MgB 2 apart from 
other diborides. 

A qualitative explanation of EFG behavior may be given based on the anisotropy of B 2p 
partial occupancies, An p =p z — {p x +p y )/2, since V zz ~ < 1/r 3 > An p and one may consider 
the boron p < 1/r 3 > expectation value to be constant for all diborides discussed. Thus, the 
variation of the electronic EFG's is determined by the interplay of p z and p x , p y occupations. 
As seen from the partial DOS (PDOS) obtained by means of a Mulliken population analysis 
, the p XtV orbitals are more occupied than are p z orbitals (Fig. 2) and V zz is negative for all 
diborides considered. For MgB 2 and BeB 2 , the high p x>y peaks lying below -2 eV lead to 
large negative An p values and, therefore, to large boron EFG's. The small increase of p z 
occupancy explains the EFG lowering for MgB 2 (and A1B 2 ) as compared to BeB 2 . Thus, 
the weaker M-B bonds for s-, p-diborides correspond to larger boron EFG's. 

For rf-diborides, the p z PDOS is more localized due to strong covalent M 3<i-B 2p bonding 
and the high peak at 3-5 eV below E^ decreases An p and V z l z compared with s-, p-diborides. 
Among the 3<i-diborides, the EFG is smallest for TiB 2 , which has the strongest p-d hy- 
bridization. Weaker p-d hybridization (less intense p^ peak) for ScB 2 and CrB 2 results in 
a larger anisotropy, An P) and the boron EFG's are larger for these compounds than for 
TiB 2 . Since Ad and 5d states are less localized than are 3d states, the corresponding B p z 
PDOS are broadened for MoB 2 and especially for TaB 2 (Fig.2) (the strong hybridization of 
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B 2p and Ta 5d states was shown also in Ref£5), that leads to the small EFG's. While a 
Mulliken analysis is not an accurate approach for the calculation of orbital charges, espe- 
cially for delocalized p-orbitals, it still allows one to describe general trends in EFG's and to 
correlate them with peculiarities of the electronic structure. Thus, we conclude that M-Bp 
hybridization is the main factor controlling the boron EFG variation. 

The effect of hydrostatic pressure on the EFG at the B site in MgB2 was investigated for 
5 and 10 GPa with lattice parameters taken from the extrapolation formulaS a = ao(l — 
0.00187P) and c = cq(1 — 0.00307P). We found a very slow increase of EFG with pressure 
(Table 1), that also demonstrates that the boron EFG's in diborides do not have a strong 
dependence on the interatomic distances. As the EFG is a very sensitive characteristic, no 
large changes are expected in the anisotropy of the B charge distribution under pressure. 

As expected, the boron p bands widen under pressure (Fig. la). One can see that the 
band shifts relative to E^ are different for different directions of BZ; the main changes in the 
occupied p x<y and p z bands occur in the low energy range at the M, K and A points. These 
bands move down with pressure relative to Ep along T-M-K-r and A-L and the overall shift 
of the PDOS to lower energies leads to the loss of these states in the energy range from E^ to 
-2 eV, as stated in Ref.El. The decrease of p X)V and p z PDOS near E^ is partly compensated 
by its increase at lower energies, and as a result, the changes in the partial p occupations are 
small. We found the increase of p occupancies with pressure to be anisotropic - the larger 
growth of p x ^ y occupancy compared to p z giving an An p increase by 0.02 for 10 GPa. The 
accurate calculation gives smaller EFG than it follows from An p . As seen from the Table, 
the increase of B p electron anisotropy is partly compensated by the increase of the core 
charge contribution, and as a result the boron EFG in MgB2 is weakly dependent on applied 
pressure. Thus, we conclude, that the charge distribution at B site shows more isotropic 
change under pressure than do the B p valence electrons, due to the compensating behavior 
of electron and core systems. NMR measurements under pressure would be important to 
confirm this conclusion. 

The PDOS changes near Ep under a pressure of 10 GPa is shown in Fig. 3. The hole 
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concentration in the bonding a bands decreases by 0.005 (within an energy interval up to 
0.8 eV). The changes in p xy PDOS is due to the behavior of these bands with pressure along 
T-A: the bands move down at V and up at A, resulting, respectively, in the loss of holes in 
the energy range up to 0.5 eV and their increase for higher energies. Although the hole p z 
PDOS change is negative in the energy interval up to 0.2 eV, its concentration is almost 
constant under pressure for energies up to 0.8 eV. The loss of both p XtV and p z states near 
Ep results in a small decrease of N(E^) by 0.02 states/eV. Thus, we showed that the carrier 
concentration decreases and so the observed decrease of the resistance with pressure is likely 
to be connected with better coupling between the sintered grains, as suggested in Ref.0. 

Finally, we estimated the pressure dependence of the Hopfield parameter, t], which is an 
electronic part of the electron-phonon coupling A = r]/M< uj 2 >, where r]= N(E^)< I 2 >. 
The calculation of the averaged electron- ion matrix element squared, < I 2 >, performed 
within the rigid muffin tin approximation!! gave a faster increase of < I 2 > with pressure 
than the N(E^) decrease. As a result, despite a small decrease of N(E^), (dN(E F )/dP = 
-0.51%/GPa), the Hopfield parameter increases with pressure as dr]/dP = + 0.55%/Gpa. 
Hence the decrease of N(E^) cannot be considered as the reason for the T c reduction, which 
is knownil to behave as dT c /dP = - 1.6 K/GPa. Thus, according to the McMillan formula, 
the main reason for the reduction of T c under pressure is the strong pressure dependence of 
phonon frequencies, which is sufficient to compensate for the electronic effects. 

Work at Northwestern University supported by the U.S. Department of Energy (grant 
No. DE-F602-88ER45372) 
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FIGURES 

Fig. 1. Band structures of (a) MgB 2 (dot lines for zero pressure, solid lines for 10 GPa), (b) 
BeB 2 , (c) A1B 2 , (d) TiB 2 , (e) ScB 2 , and (f) TaB 2 . 

Fig. 2. Boron partial densities of p X)V (dash lines) and p z (dot lines) states (PDOS) and their 
anisotropy (solid lines) for (a) MgB 2 , (b) BeB 2 , (c) TiB 2 , (d) TiB 2 . The Fermi level corresponds 
to the zero energy. 

Fig. 3. The change in the boron p XjV (solid line) and p z (dash line) PDOS in Mgi? 2 under 
pressure 10 GPa. The Fermi level corresponds to the zero energy. 
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TABLES 
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